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Abstract 

The effect of competing nonlinearity on beam dynamics in parity-time (VT) symmetric potentials is investigated. By 
using numerical methods, the existence of gap solitons is demonstrated in the first Bragg gap of optical VT symmetric 
lattices with competing nonlinearity. Meanwhile, the stability of such solitons is analyzed through introducing a small 
perturbation to the solitary solutions. The abrupt annihilation of the solitons during propagation demonstrates the 
Bragg gap solitons in VT symmetric potentials are not stable. In comparison with the on-site gap solitons, the off-site 
gap solitons exhibit more robust properties during propagation. 
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1. Introduction 

In quantum mechanics, all physical observables correspond to the eigenvalues of operators demand that the eigen- 
values should be real and thus must be Hermitian. Yet in recent years a new concept has been proposed in an attempt 
to extend the framework of quantum mechanics into the complex domain. It is found that it is in fact possible even 
for non-Hermitian Hamiltonians to exhibit entirely real eigenvalue spectra as long as they respect parity-time (VT) 
symmetry Jll 0, 01 . VT symmetry means that the eigenfunctions of a Hamiltonian are at the same time the eigen- 
functions of the PT operator, that is HPT = PTH. Generally, the action of the parity operator P is defined by the 
relations p — > —p, x — > -jc, whereas that of the time operator T by p — > -p, x — » x, i — > -/, where p and x represent 
the momentum and position operators, respectively. From this point of view, it is easy to find that a 'PT-symmetric 
Hamiltonian requires p 2 /2 + V*(-x) = p 2 /2 + V(x), which indicates that the real part of the complex potential should 
be an even function of position and the imaginary part should be an odd one. It is noteworthy to stress that this 
condition is just necessary but not sufficient. 

To date, spatial solitons (localized bound states that can maintain their shapes during propagation in a bulk 
mediafl [1] or waveguide fl) in periodic optical lattices with VT symmetry are quite involved^ H H E [H. In 



this article, we investigate the gap solitons in VT symmetric lattices with competing nonlinearity for the first time. 
The competing nonlinearity adopted in our model is the so-called cubic-quintic (CQ) nonlinearity. CQ nonlinearity 
contains two parts (proportionally to the beam intensity and the intensity square, respectively) with different signs, 
that the nonlinearity induced by the beam is greatly affected by the intensity, i.e., with different functionalities with 
respect to the intensity, the nonlinearity may change from self-focusing to self-defocusing, or from self-defocusing 
to self-focusing along the transverse profile. Even though it is reported that stable (gap) solitons are demonstrated 
in VT symmetric lattices with Kerr nonlinearity[7j, in the models combined periodic potentials 1 12] or the Bragg 
coupled-mode structure lfl3ll and the CQ nonlinearity, and in a complex Ginzburg-Landau system 1 14ll . the results we 
obtained here are quite different. The article is organized as follows: In Section 2, we briefly introduce the general 
model equation for beam propagation in VT symmetric lattices with competing nonlinearity. In Section 3, firstly, we 
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discuss the Bloch band structures of the complex potential and give the diagrams of the structures. And secondly, we 
investigate the on-site solitary solutions as well as the off-site solitary solutions in the band gap. Last but not least, in 
order to do stability analysis of the solitary solutions, we use beam propagation method to investigate the propagation 
properties of the solitary solutions with small perturbations. In Section 4 we conclude the article. 
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Figure 1: (a)-(c) Band structures corresponding to Wq = 0.2, 0.6 and 0.8, respectively. The curves are the Bloch bands and the yellow regions are 
the band gaps, (d) real parts of the eigenvalues as a function of Wo. In part I, all the eigenvalues are purely real; in part II, the eigenvalues in the 
cone-like regions are real and those in other places are complex; in part III, the eigenvalues are completely complex. The colored regions are the 
Bloch bands and the blank regions are the band gaps. 



2. Theoretical model 

We begin our analysis by considering optical wave propagation in a competing nonlinear PT symmetric potential, 
which is governed by the traditional normalized nonlinear Schrodinger equation as used in previous literaturesj7|] 

[ Tz + S + A ° [y(x) + iW(x)] ^ + ( 7 - /2 ) ^ = °' CD 

where / = \tf/\ 2 is the beam intensity, I - I 2 correspond to a competing CQ optical nonlinearitv lll5ll . Ao[V(x) + iW(x)] 
is the so called complex potential, and in our simulation we take Ao = 5. In spatial domain, the transverse coordinate 
x and the longitudinal coordinate z are scaled to the input beam width Xo and the diffraction length L&s = nokox^, 
respectively, where no is the background refractive index and Icq = 7.njA.Q. According to the necessary condition for a 
VT symmetric potential mentioned above, the real and the imaginary parts of the complexpotential should satisfy the 
relations V(-x) = V(x), W(—x) — -W(x), respectively. Similar to the previous literatureslf7l. IM ITol [Till . we consider a 
relatively simple case: 

V(x) = cos 2 (x), W(x) = Wo sin(2x), (2) 

where Wo determines the amplitude of the imaginary part and n is the period of the potential. Generally, the solution to 
Eq.([TJ has the form tfr(x) = <p{x)e il3z where <p{x) is the nonlinear eigenmode and /3 is the corresponding real propagation 
constant. Plug the solution into Eq.([TJ, we obtain 

ft + Aq[V ( x ) + iW(x)]4> + \4>\ 2 <f, - \cf>\ 4 <(> = y80. (3) 
2 



In light of the fact that <p is a complex localized wavefunction, <f> can be written as cf> — 0r + ifa, where <p R and <pi 
represent the real part and imaginary part, respectively. Substituting <p R + i<p\ for <p in Eq.Q, we end up with 
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Figure 2: (a) Power curves of solitons (solid curve for the on-site solitons and dashed curve for the off-site ones) in the Bloch gap. (b)-(c) VT 
on-site soliton field profiles for (3=1 and /? = 2, respectively, (d)-(e) PT off-site soliton field profiles for p = 1 and p = 2, respectively. The real 
parts are plotted in blue curves and imaginary parts in red curves. 



3. Bragg gap solitons 

Before solving for the localized solutions, we first analyze the linear properties of such a periodic complex poten- 
tial by omitting the nonlinear term in Eq.©. As the Floquet-Bloch theorem demonstrated, the eigenfunctions can be 
written in the following way 

<p(x) = wic(x) exp(ikx), Wk(x) = Wk(x + n), (5) 
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where k is the Bloch wave number. In this article we just consider the case that k lies in the region -1 < k < 1, i.e., the 
first Brillouin zone. Thus, we can calculate the Bloch band structures corresponding to the complex potential by using 
the plane wave expansion method (PWEI Blql . As pointed in the previous literatures [7, 8, 10], there is a critical value 
Wj? = 0.5, below which the eigenvalues are purely real. Above the threshold the eigenvalues are partially complex. 
Further increasing Wo to a certain value (the value is changing with Aq, and for Aq = 5 the value is about 0.78), bands 
will overlap each other, that indicates all the eigenvalues are complex. In FigQIa)-(c), we exhibit the three typical 
band structures for Wo = 0.2, 0.6, and 0.8. The yellow regions are the band gaps. In FigQJd), the real parts of the 
eigenvalues changing with Wo is depicted by the colored regions. The blue, red and green areas are the first, second 
and third Bloch band regions, respectively. The pink area represents the region where the first and second bands 
overlap each other. According to the critical points of Wo, the bifurcation in Fig[TJd) is divided into three parts labeled 
I (0 < Wo < 0.5), II (0.5 < Wo < 0.78) and III (0.78 < W < 1) by two dashed lines. 

Similar to the localized solutions in real potentials, in complex potentials the localized solutions can be also 
divided into two categories: the on-site lattice solitons and the off-site lattice solitons lfl7ll . And for Wo < 0.5, we solve 
the coupled Eqs.(@]i for the family of on-site as well as off-site localized solutions with real eigenvalues located in the 
Bragg gap by using the relaxation method. In Fig|2fa), we exhibit the power of the on-site (solid curve) and off-site 
(dashed curve) solitary solutions in Bragg gap (0.54 < /3 < 2.78) for Wo = 0.2. From FigJSJa), firstly we can see that 
the power of the off-site solitons is bigger than that of the on-site ones, and secondly we can see that the power for 
both two type solitons has a minimum value at /3 very close to the right boundary of the band gap. Corresponding to 
the red dots in Fig|2la), we display the localized solitary solutions in Figj2^b)-(e). Comparing Fig|2lb) with Figj2jc) 
and Fig|2d) with Fig|2e), and considering the power curves shown in Fig|2a), we can conclude that (i) the more 
localized the solitary solution is, the lower the corresponding power is; (ii) very close to the boundaries of the Bragg 
gap, it is really hard to observe a solitary solution, for the localization becomes worse and worse. 




Figure 3: (a) The intensity of the on-site gap soliton (blue curve) and the index change induced by the soliton (red dashed curve); (b) Propagation 
of the soliton shown in (a) with a small perturbation; (c) The intensity of the off-site gap soliton (blue curve) and the index change induced by the 
soliton (red dashed curve); (d) Propagation of the soliton shown in (c) with a small perturbation. The parameters are Wo = 0.2 and/? = 1 for both 
solitons. 

As did in the previous stability analysis of the solitons in PT symmetric potentials discussed in this 

article is also very important. In Figf3] we show the propagation dynamics for the two soliton solutions displayed in 
Fig|2fb) and (d) with small perturbations both on the amplitude and phase. It is clear to see that the beams annihilate 
abruptly if the propagation distance exceeds a certain value, which is similar to the properties of a "soleakon" reported 
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in 2009 [20], but they do not share the same physical mechanism. One main difference is that the input used by a 
soleakon is leaky mode which is unbound state, while the input used in this article is really bound state. The reason 
why the solitons are not stable lies in the fact that defocusing nonlinearity does not support stable bright solitons. 
As shown in FigfJJa) and (d), the index changes (red dashed curves) induced by the humps of the gap solitons (blue 
curves) are negative, that means the nonlinearity is defocusing. Anyway, in a long range (e.g., < z < 150 for the 
case shown in FigJ3b) the solitons are relatively stable. In a physical point of view, if we take a soliton beam with 
Aq =1.0 fim, rtQ = 1.5, and xq = 10 //m as the incidence, it can stably propagate about IkoHQX^z ~ 28 cm with a 
diffraction length Lditr ~ 1 mm. In comparison with the on-site gap solitons as shown in Fig|3jb), the off-site gap 
solitons as shown in Figj3jd) are much robust. 
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Figure 4: (a) Bragg gap diagram changing with Aq (the parameter Wo = 0.2), in which the semi-infinity gap and Bragg gap regions are marked, 
(b-c) The solitary solution (left panel) and its propagation (right panel) corresponds to the blue dot. (d-e) The solitary solution (left panel) and its 
propagation (right panel) corresponds to the green dot. 

Naturally, people may ask can the Bragg gap solitons in the complex potentials with higher amplitude Aq be stable 
during propagation? To answer this question, we first display the Bragg gap diagram changing with Aq as shown in 
FigHJa), in which the two red dots correspond to the unstable solitons we discussed in Figf3](Ao = 5). And then, we 
discuss another two cases at Ao = 10 (blue dot) and Aq = 20 (green dot), as shown by Figsj4jb-c) and FigsHJd-e), 
respectively. According to the propagations, we can conclude that the answer to the question is negative, for after a 
relatively stable and breath-like propagation, the Bragg gap solitons annihilate eventually. 



5 



4. Conclusion 

In conclusion, the existence of Bragg gap solitons in periodic VT symmetric potentials with CQ nonlinearity 
is demonstrated numerically. Even though the gap solitons we found are not stable, for they will annihilate during 
propagation because of defocusing nonlinearity they induced, they can still stably propagate over a long distance. 
Contrast of the two propagations of on-site solitons and off-site solitons, we find that the latter ones are more robust. 
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